Energy relaxation of a superconducting charge qubit via Andreev processes 
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We study fundamental limitations on the energy relaxation rate of a superconducting charge 
qubit with a large-gap Cooper-pair box, Af, > A r . At a sufficiently large mismatch between the gap 
energies in the box A& and in the reservoir A r , "quasiparticle poisoning" becomes ineffective even in 
the presence of nonequilibrium quasiparticles in the reservoir. The qubit relaxation still may occur 
due to higher-order (Andreev) processes. In this paper we evaluate the qubit energy relaxation rate 
Tj -1 due to Andreev processes. 



I. INTRODUCTION 

A large number of recent experimental stud- 
ie ol'2i3i4i5iEiL9i? indicates the presence of quasiparticles 
in superconducting single-charge devices at low tem- 
peratures. The operation of these devices, of which 
the best known is Cooper-pair box qubit, requires 2e- 
periodic dependence of the charge of the box on its 
gate voltage, and thus, an introduction of an unpaired 
electron(quasiparticle) in the Cooper-pair box (CPB) 
is a significant problem. The superconducting charge 
qubit operates at the degeneracy point for Cooper-pairs, 
N g = 1, with N g being the dimensionless gate voltage. 
For equal gap energies in the Cooper-pair box and reser- 
voir, Ab = A r , the states of the qubit at N g = 1 are 
unstable with respect to quasiparticle tunneling to the 
box. The quasiparticle changes the charge state of CPB 
from even to odd, and lowers the charging energy. This 
phenomenon, commonly referred to as "quasiparticle poi- 
soning" , is well-known from the studies of the charge par- 
ity effect in superconductors, see, for example, Matveev 
et. alJ£ and references therein. "Quasiparticle poison- 
ing" can degrade the performance of the charge qubit 
in two ways. First, it causes the operating point of the 
qubit to shift stochastically on the time scale comparable 
with the measurement time^. Second, it contributes to 
the decoherenceii. One of the approaches to improve the 
performance of charge qubits is to use superconducting 
gap engineering. In most single-charge superconducting 
devices "quasiparticle poisoning" can be suppressed even 
in the presence of nonequilibrium quasiparticles in the 
reservoir by engineering a large mismatch between Af, 
and A r . Gap energies in superconductors can be modi- 
fied by oxygen doping^, applying a magnetic field^, and 
adjusting layer thickness^. In this paper we study the 
fundamental limitations on the energy relaxation time in 
a charge qubit with a large gap in the box, A& > A r . 

For equal gap energies in the box and reservoir, A& = 
A,., the energy relaxation rate due to "quasiparticle poi- 
soning"— is 



T\ hi>p \/ E, 



(1) 



tion and density of states at the Fermi level, respectively. 
The relaxation rate 1/Ti in Eq. |T]) was derived under the 
assumption that an unpaired electron tunnels from the 
reservoir to the box to minimize the energy of the system. 
Indeed, for A;, = A r , the odd-charge state of the CPB 
has lower energy at N g = 1 due to the Coulomb block- 
ade effect. By properly engineering superconducting gap 
energies (i.e. inducing large gap mismatch, A& > A r ), 
one can substantially reduce quasiparticle tunneling rate 
to the Cooper-pair box. Suppose initially the qubit is in 
the excited state with energy and the quasiparti- 

cle is in the reservoir with energy E p . Upon quasipar- 
ticle tunneling to the box, the minimum energy of the 
final state is Ef m = A5 + En+i with En+i being the 
energy of the CPB in the odd-charge state. Therefore, 
the threshold energy for a quasiparticle to tunnel to the 
box is £™ in = A & + E N+1 - E\ + ), see also Fig. Q] If 
Ep Un — A r > E j ^> T, only exponentially small fraction 
of quasiparticles are able to tunnel into the island. (Note 
that the energy difference between excited and ground 
state of a charge qubit is E,, while the energy of the 
qubit in the excited state is E\+\ — E c + E ; /2. Here 
E c , E } and T are the charging energy of the CPB, the 
Josephson energy associated with the tunnel junction, 
and the temperature, respectively.) Thus, the contribu- 
tion to the qubit relaxation rate T^ 1 from the processes 
involving real quasiparticle tunneling to the island be- 
comes 



1 #T n qp 

— oc — — — exp 



-A r -E c -Ej2 



(2) 



with n qp , g T and vp being the density of quasiparticles 
in the reservoir, dimensionless conductance of the junc- 



and is much smaller than the one of Eq. (fT|). (To obtain 
Eq. we used the fact that -Ejv+i = at N g = 1.) 
However, there is also a mechanism of energy relaxation 
originating from the higher order tunneling processes 
(Andreev reflection) . The contribution of these processes 
to the qubit relaxation is activationless, and can be much 
larger than the one of Eq. ([2]). In the rest of the paper we 
study qubit energy relaxation due to Andreev processes 
in detail. 
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FIG. 1: (color online). The spectrum of the Cooper-pair box 
as a function of the dimensionless gate voltage for a large-gap 
mismatch, A;, > A r . The solid and dashed lines correspond 
to an even- and odd-charge state of the box, respectively. 



II. THEORETICAL MODEL 

Dynamics of the Cooper-pair box coupled to the su- 
perconducting reservoir through the tunnel junction is 
described by the Hamiltonian 



H = H r 



(3) 



Here H^ cs and H^ cs are BCS Hamiltonians for the box 
and reservoir; H c = E c (Q/e — N g ) 2 with N g and Q be- 
ing the dimensionless gate voltage and the charge of the 
CPB, respectively. We consider the following energy scale 
hierarchy: Ab > A r ,E c > E , ^> T . In order to distin- 
guish between Cooper-pair and quasiparticle tunneling, 
we present the Hamiltonian ([3]) in the fori 



,12 



H = H + V, and V = H^-H 



(4) 



Here H = H, 



Hbcs 



■m 



BCS 



-H , , and H , is the Hamil- 



tonian describing Josephson tunneling 



H, 



\N) (N\ H T 



1 



E-H 



H T \N+2) (N+2\+R.c. 



The matrix element (N\ H T ^ r H T \N+2) is propor- 
tional to the Josephson energy E r . The perturbation 
Hamiltonian V defined in Eq. (0| is suitable for calcula- 
tion of the quasiparticle tunneling rate. The tunneling 
Hamiltonian for homogeneous insulating barrier is 

Ht = J2 I rfxdx ' (7 1 (x,x')^(x)* ff (x')+H.c.) , (5) 

a J 

where x and x' denote the coordinates in the CPB and 
reservoir, respectively, and T(x, x'), in the limit of a bar- 
rier with low transparency, is defined as 



T(x,x') 




6 2 (r-r / )S(z)S(z') 



,.dd 



dz dz' 



(6) 



Here T is the transmission coefficient of the barrier, r and 
z are the coordinates in the plane of the tunnel junc- 
tion and perpendicular to it, respectively. The Hamil- 
tonian ([5]) along with the above definition of T(x, x') 
properly takes into account the fact that in the tunncl- 
Hamiltonian approximation the wavefunctions turn to 
zero at the surface of the junctio n 13 ' 14 . In terms of 
the transmission coefficient T, the dimensionless con- 
ductance of the tunnel junction g T can be defined as 
g T = TSjkp/Air = ^TN ch , where Sj is the area of the 
junction, and is the number of transverse channels 
in the junction. 

The energy relaxation rate of the qubit due to higher- 
order processes is given by 

r^ = yE 2 \ A *p\*8( E P>- E p- E j)fAEp)(l-fA E P>))- 



p,p' 



(7) 



Here f F (E p ) is the Fermi distribution function with 
E p = ,/e^ + A 2 . being the energy of a quasiparticle in 



the reservoir. The amplitude A p i p is given by the second 
order perturbation theory in V, 



Ap>p = (-, Ep>i \ V 



1 



Ei — H 



V\+,E p] ) 



(8) 



At E c 3> Ej and N g = 1, the eigenstates of the qubit 
are given by the symmetric and antisymmetric super- 
position of two charge states, i.e. | — ) = \ N )+\^ +2 ) 



and 



\N)-\N+2) 
V2 



with the corresponding eigenval- 
ues E\±\ — E c ± Ej/2. In the initial moment of time the 
qubit is prepared in the excited state and the quasipar- 
ticle is in the reservoir, i.e \+,E p f) = |+) ® \E p ^). The 
energy of the initial state is Ei — E p + Ei + \ . The denom- 
inator in the amplitude ((5|) corresponds to the formation 
of the virtual intermediate state when the quasiparticle 
has tunnelled to the island from the reservoir. Since a 
quasiparticle is a superposition of a quasi-electron and 
quasi-hole, the contributions to A p i p come from two in- 
terfering paths: 



-2,E P ^\V 



1 



Ei — H 



V\N,E pr ) 



2 W^WI 



V- 



1 



Ei — H 



V\N+2,E, 



(9) 



To calculate the amplitude A p ' p , we use particle- 
conserving Bogoliubov transformatio n 15 ' 16 ' 17 : 



7 t 

'7 la 



7n<r 



J dx [C/„(x)*t ( x ) _ CT V„(x)*_ a (x)iit] 

f7„(x)* ff (x) - CT^ n (x)*L ff (x)J2l (10) 



= / dx 



The operators W and R transform a given state in an N- 
particle system into the corresponding state in the N + 2 
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and N—2 particle system, respectively, leaving the quasi- 
particle distribution unchanged, i.e. RJ \N) = \N + 2). 



Thus, quasiparticle operators 



7 t 

/ ncr 



and 



In 



defined in 



Eq. (JTU]) do conserve particle number—. The transfor- 
mation coefficients U n (x) and V^(x) are given by the 
solution of Bogoliubov-de Gennes equation. For spa- 
tially homogenous superconducting gap A, the functions 
L/„(x) and V n (jx.) can be written as Z7 n (x) = it„0 n (x) and 
V n (yi) = v ra </>n(x). The coherence factors u n and v n are 
given by 



values and eigenfunctions of the single-particle Hamilto- 
nian, which may include random potential V(x), e.g., due 
to impurities. The single-particle energies e„ and wave- 
functions 0„(x) are defined by the following Shrodinger 
equation: 



2m 



V 2 + V(x) 



</>n(x) 



»(x) 



I" 



) and vl 



2 I ' E„ 



Here E n = \/e* l + A 2 ; e n and </>„(x) are exact eigen- 



In the presence of time-reversal symmetry u n , v n and 
</> ra (x) can be taken to be real. Then with the help of 
Eq. (fTU|) , we obtain the amplitude of the process A p i p : 



\ J dxidxidx^TtxiXmx^) [U P ,(A)V P {A)-U P {A)V P >{A)\ E eI+6E+-eI 



(11) 




b) 



Reservoir 



Box 




Reservoir 



FIG. 2: The diagrams corresponding to the interference of 
electron trajectories in the box (a) and reservoir (6). The 
contribution of the diagrams with interference in both elec- 
trodes (not shown) is much smaller than the one of the above 
diagrams^. 



where SE + = Ei + \ — E^+x = E c + E ,/2. The minus 
sign in the parenthesis here reflects the destructive inter- 
ference between quasi-electron and quasi-hole contribu- 
tions, see also Eq. ©. 



III. DISORDER AVERAGING 

It is well-known that Andreev conductance is sensitive 
to disorder, see, for example, Refs. 1_9,20]. Similarly, the 
rate Ta is affected by electron backscattering to the tun- 
nel junction, see Fig. [2j If a quasiparticle bounces off 
the walls of the box or impurities many times, it is rea- 
sonable to expect the chaotization of its motion. Thus, 
one is prompted to consider ensemble-averaged quantities 
rather than their particular realization. Using Eqs. (O 
and pip , we obtain 



J 



(T A ) = / II dx i rfx i r ( x l I x i) T ( x 2,x' 2 )T(x3,X3)T(x 4 ,X4) (Up'Vp^ p '(-X.' 1 )(j)p(x' 2 )- 

p,p' 2—1. .4 



-U p Wp'0p(x / 1 )(/) p / (x' 2 )) 



X (Up,V p <pp,{x 3 )(t>p(x4)-U p Vp,<f>p(x 3 )(l> p ,{X4)) 2^ - +SE+ _ Ek 2^ ■ 



E p + 5E + — Ek> 



x 5(E p , -E p - Ej)f F (Ep)(l - f F (E p ,))). 



(12) 



Here the brackets (...) denote averaging independently in the CPB, one has to calculate the following correlation 
over different realizations of the random potential in the 
box and reservoir. In order to average over the disorder 



4 



function: 



7= 



E 



U k Vk 4>k (Xi ) 0fc (X 2 ) U k > V k > 4>k> (x 3 ) </>fc' (x 4 ) 



\k,k' 



E P + SE+ — E k 



E p +5E+-E k , 



■I 



(K (l (x 1 ,x 2 )K S2 (x 3 ,x 4 )) 



4E(£ 1 )E(ta)[E p +6E + -E(t 1 )][E p +5E+-E(t a ) 



(13) 



where JtTg(x 1 , x 2 ) = £) fc <fc(xi)0fc(x 2 )5(efe - £), and 
A(£) = ^£ 2 + A 2 . The correlat ion function 
(if{ 1 (xi,x 2 )A^ 2 (x 3 ,x 4 )) consists of reducible and irre- 
ducible parts, 

(if Sl (xi , x 2 ) Afr 2 (x 3 , x 4 )) = (K £l (xi , x 2 )) (iQ 2 (x 3 , x 4 )) + 
+ (if a (xi,x 2 )^ 2 (x 3 ,x 4 )) ir . (14) 

The reducible part can be easily calculated by relating 
(-?Q(xi,x 2 )) to the ensemble-averaged Green function: 
(JQ(xi,x 2 )) = -ilm(G| f (xi,x 2 )) = v F fi 2 . (Upon av- 
eraging over disorder, one can neglect the energy depen- 
dence of the density of states here, i.e. (vf(Q) = Vf- 
The function /12 is given by /12 = (e lk ( xl ~ X2 ))FS with 
(...)fs being the average over electron momentum on 
the Fermi surface. For 3D system the function /12 is 
equal to /12 = S "fcpfxi-X2 \ •) ^he irreducible part 
(K^ (xi, x 2 )A~£ 2 (x 3 , x 4 ))i r can be expressed in terms of 
the classical diffusion propagators - diffusons and Cooper- 
ons, see, for example, Aleiner et. cd. [2l|. In the ab- 
sence of magnetic field, diffusons and Cooperons coin- 
cide, 7^(xi,x 2 ) = 7^(xi,x 2 ) = 7 3 5(xi,x 2 ), and the 
irreducible part of the correlation function (| 14[) reads 

(A Cl (xi,x 2 )A^ 2 (x 3 ,x 4 )) ir = (15) 

= ~ RC [/l4/23^|C 2 ^ 1 |(xi,X 3 ) + / 1 3/24'P| ?2 ^ 1 |(xi,X4)] . 

The spectral expansion of V u (x! , x 2 ) for the diffusive sys- 
tem is 



7 ? w (xi,x 2 ) = 



/,t(x 1 )/»(x 2 ) ^ 
-iuj + 7„ 



(16) 



Here "f n and / n (x) are the corresponding eigenvalues and 
eigenf unctions of the diffusion equation, — 7>V 2 /„(x) = 
7n/n(x), satisfying von Neumann boundary conditions in 
the box. 

Equation (fT5|) can be simplified in the case of large 
Thouless energy, i.e. E T > A b , A r , E c , E s . (Here E T — 
K/td with td ~ Sb/D being the time to diffuse through 
the box, and Sb being the area of the island, see Fig. [3]) 
This condition is fulfilled for a small aluminum island 22 
with St <C 1/xm 2 and mean free path I > 25nm 2 ^, when 
the time spent by the virtual quasiparticle in the box, 
t ~ Kj (Ab — A r — 8E+), is much longer than the classical 



diffusion time tet^. In this case the irreducible part in 
Eq. (fT4"|) is given by the universal limit, 



( K Si ( x i , x 2 ) A ?2 (x 3 , x 4 )) ir 



(17) 



V F 



= TT<*(£l - 6)(fu/23 + /l3/24)- 
Vb 

Here Vb is the volume of the box. Upon substituting 
Eqs. (HU) and ([TTJ) into Eq. (TJ3]) and evaluating the inte- 
grals over energies £1 and £2, we obtain 



I = 414/12/34A1 



E P + 5E+ 



^2A~~ ^ 14 ^ 23 + ^ 13 ^ 24 ) ^ 2 



E P +5E+ 



(18) 



where Sb = 1/vpVb is mean level spacing in the box. The 
functions Li(y) and L<z(y) are defined as 



Lx{y) = 



1 



■ arctan 



1+2/ 



1 1 

Vx 2 -1 x(a;-y) 2 



(19) 



The expressions above are valid for y < 1. The function 
7/2 (2/ ) has the following asymptotes 



Ml/)' 



- + -v 

4 ' 3 y > 



I/<1, 

r, 1-»<1. 



(20) 



2 v / 2(l-a) 3 / 2 ' 

After substituting Eq. (fTg)) into Eq. (TT2]) and averaging 



CPB 



Reservoir 




FIG. 3: The layout of the Cooper-pair box qubit considered 
in the text. 

over disorder in the reservoir, we obtain the following 
expression for (r^): 
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{Ta) ^ I ^i^W 2 )-^i)-^)/fWi)](1-/f[B(^)])/ II rfx 4 ^T(x ll x' 1 )T(x 2 ,x^)T(x3,x^)T(x 4 ,x^) 



x 4/12/34X1 



A;, 



2A, 



(/14/23 + /13/24) £2 



i=1..4 



A,, 



(21) 



Here = \/£' 2 + A 2 . The correlation function 

in the reservoir (Kg (x'i, xfg^Kg (x'2, x'4)) follows from 
Eqs. (fl4|) and (fPoj). Using Eq. (O and evaluating the 
spatial integrals over the area of the junction as well as 
the integrals over energies £[, and £' 2 , we finally obtain 
the answer for (Fa): 



(r A ) = rx + r 2 

with Ti and T 2 being defined as 



E, 



2tt 3Ci g 2 T 

ft (4^ 2 ) 2 7V ch \/2A r + J B J ^ F 



A r +(5£' + 



A?, 



and 



2tt 



h 8(4tt 2 ) 2 A b V 2A r +J5 J v F 



-SE 4 



(22) 



,(23) 



.(24) 



Here C\ is a numerical constant of the order of one: 



Ci 



1 



dyi dy 2 dy 3 dy 4 P 12 P 13 P 24 P34 



with y being a dimensionless coordinate in 
the plane of a tunnel junction, and Pyi = 

sin(| y i- y2 |)-|yi-yj,|coB(| yi - y2 |) Th functions L 

|yi-y2| 3 1 
and L 2 are defined in Eq. (|19p . and their dependence on 
the ratio (A r + <5£'+)/Ab is shown in Fig. 0J The rate 
Y\ describes the contribution from the reducible terms, 
see Eq. (fT4"]) . and is similar to the ballistic case when 
electron scattering from the impurities or boundaries is 
negligible. The other term, T 2 , reflects the enhancement 
of (r^) in the diffusive limit due to the quantum 
interference of quasiparticle return trajectories^, and 
originates from the irreducible contributions, see Fig. [2] 
In the case of N^Sb/Ab ^> 1, the contribution of this 
interference term becomes dominant, T 2 Ti. The 
contribution of the interference in the reservoir to the 
rate r 2 , see Fig. [5b, is geometry dependent. For a 
typical charge qubit with the small junction connected 
to a large electrode, backscattering of electrons to the 
junction from the reservoir side gives much smaller 
contribution to T 2 than the similar one for the box side 
of the junction. In particular, for the layout of the qubit 
shown in Fig. [3l the contribution of the interference 
in the reservoir to T 2 is smaller than the one in the 



box by a factor ^ At m 



hD 
A r Sj 



Li(a ) 
£-2(00) 



< 1. [Here 



L{a ) 




FIG. 4: The dependence of the functions Li(ao) and £2(00) 
(normalized by Li(0) and 1/2(0), respectively) on the dimen- 
sionless parameter ao = (A r + 5E+)/At- The solid and 
dashed lines correspond to L\ and L2, respectively, and re- 
flect the increase of the rates Fi and 1*2 with ao. The expres- 
sions for Li(ao) and I/2(ao) given by Eq. (|19[) are valid for 
00 < 1 - T/A b . 



ao = (A r + 8E+)/ Ab, and db( r ) is the thickness of the 
superconducting film in the box(reservoir).] Therefore, 
we neglected the terms corresponding to the interference 
in the reservoir in Eq. (|24p . 



IV. CONCLUSION 

We have studied the fundamental limitations on the 
energy relaxation time in a charge qubit with a large- 
gap Cooper-pair box, Ab > A r . For sufficiently large 
Ab, real quasiparticle transitions can be exponentially 
suppressed, and the dominant contribution to the charge 
qubit energy relaxation time T\ comes from the higher- 
order (Andreev) processes, see Eq. (|2"2"1) . For realistic ge- 
ometry of the charge qubits and the density of nonequi- 
librium quasiparticles in the reservoir n qp ~ 10 19 — 
10 18 m -3 [ll|, we estimate the Andreev relaxation rate 
to be (r^) ~ 1CT 1 - l(T 2 Hz. Thus, in the absence of 
other relaxation channels, the mismatch of gap energies 
leads to extremely long Ti-times. (For comparison, the 
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quasiparticle-induced Ti found in Ref. [ll[ for the charge 
qubit with equal gap energies was T^ 1 ~ 10 5 — 10 3 Hz.) 

The charge qubit with a large gap in the box also per- 
mits to reduce quasiparticle-induced decoherence. Since 
real quasiparticle transitions into the island are sup- 
pressed, see Eq. ([2]), the dephasing time of the qubit is 
limited by the energy relaxation processes, i.e. T2 ~ 
2/<Ia>. ' 
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